Damping mechanism for the strongly renormahzed c-axis charge transport 
in high-Tc cuprate superconductors 
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One of the peculiar aspects of the high- Tc cuprate su- 
perconductors is the large anisotropy related to their lay- 
ered structure. Although the anisotropy following from 
a conventional Fermi-liquid approach based on LDA (lo- 
cal density approximation) band structure calculations, 
is of the order 10 in most of these compounds, experi- 
mentally it is found to be much larger, as follows from 
e.g. the London penetration depth and the normal-state 
resistivity. In the clean limit the London penetration 
depth Al is given by c^/A| = Anne'^/'m* and usually 
the large anisotropy of Al is attributed to a large mass 
enhancement ||^,|^. A mass enhancement could follow 
from e.g. spin-charge separation [^JJJ^. Alternatively, 
elastic g or inelastic [Q Fermi-liquid corrections, as well 
as spin-charge separation, may cause a strongly incoher- 
ent charge transport along the c axis. No plasma pole 
is found in the normal state, whereas in the supercon- 
ducting state in some of the materials a zero-crossing of 
the real part of the dielectric function (e') is found at 
an energy two orders of magnitude below the value ex- 
pected from LDA band structure calculations. In this pa- 
per we shall argue that the discrepancy results not from 
a strong mass-enhancement, but from a strong quasi- 
particle damping for E//c, both in the normal and in 
the superconducting state. 

Single crystals of La2-:ESra;Cu04_5 with nominal Sr 
contents of 10 and 15 % were grown by the traveling- 
solvent floating-zone method Electron probe micro 
analysis showed the correct phase with no contamina- 
tion or inclusion of second phases, and a homogeneous 
distribution of strontium. The crystals have sharp su- 
perconducting transitions (35 K at optimal doping, with 
90 % of the transition within 3 K as determined from AC 
susceptibility measurements at 3 Gauss, 90 Hz). 

Neutron diffraction measurements demonstrated that 
the as-grown boules contain large single crystalline 
grains. Samples for reflectivity measurements were cut 
perpendicular to the b axis and were mechanically pol- 
ished to optical quality. The surface in the ac plane 
was approximately 4x7 mm^. X-ray Laue backscatter- 
ing showed a correct orientation with sharp, well-defined 
single diffraction spots. 

Polarized infrared reflectivity spectra were measured 
in two different laboratories between 10 and 10000 cm~^ 
at temperatures between 10 and 300 K using rapid- 
scan Fourier-transform spectrometers equipped with a 
continuous- flow cryostat. The data obtained in the two 
laboratories on nominally identical crystals are in a good 
agreement with each other, and with those of Tamasaku 
et al. 0. 

In the inset of Fig. 1 we display the reflectivity of 
a Lai gsSrg 15CUO4 single crystal with E//cm the nor- 
mal state at 33 K. The spectrum is typical for an ionic 
insulator. At low frequencies there is a Hagen-Rubens 
type of upturn, i?(t^) = 1 — 2((To/j^)^/^, characteris- 
tic of a poor conductor, where Iiea{iy) is dominated by 



cro(l + with tTo ~ 7 n~-^cm~^. A closer in- 

spection, using a Drude-Lorentz model with a Drude os- 
cillator and 4 phonons (see Fig. 1 caption), indicates a 
featureless and almost constant value of the electronic 
part of (t(i'): A good fit in the range of 20 to 650 cm~^ 
can only be obtained (see Fig. 1 inset) if 7 is taken to 
be larger than 0.5 eV. With these parameters a zero- 
crossing of e' = £00 — 2(To7(i^^ +7^)~^ occurs for v"^ < 0, 
and hence the plasmon has a diffusive pole in the normal 
state. By integrating the corresponding expression for 
Recr(j/), we can put a lower bound on Up of ^/2aoJ = 0.15 
eV. Interestingly this is much closer to estimates of the 
c-axis plasma frequency based on the single-electron ef- 
fective mass in this solid. From such an RPA-LDA anal- 
ysis for the La2Cu04 system one obtains a screened ab- 
plane plasma frequency of 2 eV, well above the onset of 
the interband electron- hole continuum at 0.5 eV, and a 
screened c-axis plasma frequency at 0.35 eV, i.e., below 
the interband electron- hole continuum It is interest- 
ing to notice that experimentally the situation vis d vis 
Landau damping is reversed: In the a6-plane dielectric 
function a clear, though strongly damped, plasma pole is 
observed near 1 eV | [lO| |, whereas a plasmon is absent in 
the c-axis dielectric function due to overdamping. Such 
loss of coherence is also born out by confinement of holons 
to the planes [|ll| , due to the fact that an interlayer hop- 
ping process requires a spinon-holon combination as an 
intermediate state. 

Let E* be the effective value of the gap, or gap- 
distribution. We will assume E* = 3.5A:bTc, which be- 
comes exact in the limit of a weak coupling s-wave super- 
conductor. In the superconducting state we should make 
a distinction between the high- and low-frequency parts 
of the spectrum, far above and below E* (w 80 cm~^ for 
Lai.85Sro.i5Cu04). The former remains relatively unaf- 
fected by the transition, whereas for the latter range of 
frequencies we expect to find temperature dependence 
in both the real and imaginary parts of the electronic 
contribution to the conductivity. For example, for an 
isotropic s-wave BCS superconductor one expects a de- 
crease of a{i') for v K E*/2 as indicated in Fig. 2. In 
addition there may be a strong decrease of 7 for T < 
as has been reported for the in-plane infrared [p^ , p^ and 
microwave response of high-Tc superconductors, and 
the c-axis infrared response of La2-xSr2.Cu04 Q. 

We analyze our data in the superconducting state using 
the two-fluid Gorter-Casimir model at frequencies up to 
100 cm~^, i.e., with the model dielectric function [nsl 



+ i 



2an 



1^(1 — iv 1^') 

The dashed curve in Fig. 1 corresponds to the fit- 
parameters used by Tamasaku et al. for 8 K 0|. Note 
that with their set of parameters in the superconduct- 
ing state nearly all oscillator strength under a collapses 
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into the 5-function at the origin. This essentially cor- 
responds to the clean limit. A much better fit is ob- 
tained if the spectra are modeled with a broad and fea- 
tureless electronic component of cr, similar to the normal 
state. Initially we varied all parameters. At all tem- 
peratures the best fits indicate that vj^ ^ 1, and only 
small and unsystematic changes of Eoo with temperature 
are obtained. Therefore we repeated the analysis with 
7 — + oo, and replacing Eqc by the sum of Coo and the 
oscillator strengths of Fig. 1. To estimate the influ- 
ence of dispersion of the optical phonons, we also did 
this analysis using the full set of phonon parameters of 
Fig. 1. The temperature dependence of the remaining 
two parameters ((i^0(T)/i^0(O))^ and cro(r)/(To(Tc), with 
1/^(0) = 180 cm"-'^ and (To(31 K) = 6.7 ri"^cm~^ are dis- 
played in Fig. 2, and the reflectivity data together with 
the fitted curves are shown for three different tempera- 
tures in Fig. 1. 

Also shown in Fig. 2 are the results of model cal- 
culations for the dirty limit of isotropic s-wave 
and d-wave superconductors, assuming standard ra- 
tios of E*/Tc for these two models, and assuming a 
temperature-independent scattering rate. The conduc- 
tivity was calculated for v = cm^^. This frequency 
is close both to BCS ratio S.^ksTc — 80 cm~^ and to 
the gap value deduced from neutron scattering [|8|. As 
a result we observe a steep drop of a for T < Tc, and the 
change of analytical behavior at the temperature where 
Eg(T) = hv. For d-wave pairing the calculated tem- 
perature dependence of cr is smooth, and quite close to 
the experimentally determined values. If we would as- 
sume that in addition there is a drop of the decay rate 
when the gap opens, the theoretical curves would have 
an even steeper temperature dependence below Tc. For 
this reason we attribute the decrease of a for T <Tc ex- 
clusively to the opening of a gap, i.e., without assuming 
a change of the electronic decay rate. For (To(r) a rea- 
sonable agreement exists with the d-wave model, whereas 
there is no good agreement with the s-wave model. For 
fs{T) = {i>^{T)/u^{Q)Y the situation is reversed: fs fol- 
lows roughly a 1 — (T/Tc)* law, which differs strongly 
from the d-wave pairing model, and lies somewhat above 
the s-wave pairing curve. For cro(r) the comparison with 
the s-wave model is slightly improved if we postulate 
the existence of constant background (possibly spurious) 
conductivity in our samples. Interestingly, the jump in 
X^(r)/x"''(r), as recently reported by Shibauchi et al. 
1^, is very well reproduced within this picture, as well as 
the difference between Xab and Ac, if we assume a clean 
limit picture for the a6-plane response. 

From a direct inspection of (7(1^) obtained from 
Kramers-Kronig analysis (Fig. 1 lower panel) we find 
no signature of a gap in the data ||l^, possibly due to 
limitations in signal-to-noise ratio and the presence of a 
tail of the optical phonon at 240 cm"-'^. As is clear from 
the calculated conductivity (the dashed curve) based on 



Mattis-Bardeen theory, even in the case of an isotropic s- 
wave superconductor a gap would be hard to distinguish. 
In the present case the difficulty arises from the low value 
of a due to a strong damping. We notice that the sit- 
uation is reversed compared to the planar response. In 
the latter case E* is situated in the Lorentzian tail of a 
narrow Drude peak, which, due to the low value of cr, 
makes a gap difficult to observe pOj . 

In Fig. 3 we classify various regimes that can occur de- 
pending on the relevant energy scales: First there is the 
pair-breaking energy scale E*. Second is the screened 

plasma frequency Vp — tat^'^Vp, (the screening may in- 
clude optical phonons depending on how far down the 
plasmon is pushed due to screening and damping effects). 
Similar parameters were also used by Uemura et al. p^ ] 
to classify clean-limit superconductors corresponding to 
their Fermi energy (oc hvp for a 2D Fermi liquid). We 
extend this classification with the the scattering rate 7 
(using 270-9 = Vp). We scale both E* and 7 to 

The region labeled 'CL' corresponds to the clean limit 
where j E*, and this region contains the 'exotic' su- 
perconductors considered by Uemura et al. p^ , in par- 
ticular the high- Tc cuprate superconductors character- 
ized by their in-plane electrodynamical properties. Most 
classical superconductors fall in the region labeled 'MB', 
where the infrared properties can be described in the con- 
text of dirty-limit Mattis-Bardeen theory. In the CL and 
MB regions a plasmon is found well above E* in the nor- 
mal as well as in the superconducting state. 

In the regions labeled 'OP' and 'RP', the plasmon is 
overdamped in the normal state, as with ^ > Up a zero 
crossing of e' no longer exists. In the superconducting 
state the 'missing area' of the conductivity due to the 
opening of the gap is transferred to the (5-function at 
J/ = 0, so that e' crosses zero at Vp^'KE*l2h'^. Un- 
less this frequency is below the gap, the plasma mode 
will be damped (OP region). Hence the requirement 
for the existence of a long-lived collective mode is, that 
E*/vp > TTVp/2^. This 'phason' mode has been described 
by Doniach and Inui | |2^ assuming Josephson coupling 
between the layers. Using Anderson's expressions |2^] 
for collective modes in the superconducting state we have 
recently shown that the plasmon gradually changes 
from a collective oscillation of quasi-particles for ly E*, 
to an oscillating current carried by the condensate for 
V <^ E*. This region is labeled 'RP' (re-entrant phason). 
Our investigation indicates that the c-axis response of 
Lai.85Sro.i5Cu04 lies in the RP region. Although the c- 
axis response data for other high- Tc materials are very 
limited, the far infrared results of Noh et al. |2^ and of 
Tajima et al. [ p7| indicate that also the c-axis response of 
Ndi.g5Ceo.i5Cu04-y and of Bi2Sr2CaCu208-i-£!; should be 
placed in the RP region. In both these materials the c- 
axis plasmon, attributed to a sphere resonance by Noh et 
al. [^, is found at very low energy of approximately 10 
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cm~-^, and in Bi2Sr2CaCu2 08+a; a very strong damping 
of c-axis electronic response is observed. 

The two regions labeled CLP ('clean local pair') and 
OLP ('overdamped local pair') correspond to clean and 
dirty limit superconductors with a pair-breaking energy 
larger than the plasma frequency. In such systems the 
phase transition is induced by thermal activation of pha- 
sons, while the pairs remain intact up to a higher de- 
pairing temperature scale. This corresponds to Bose- 
Einstein condensation of local pairs. According to the 
analysis of Tamasaku et al. the c-axis phason of 
La2-2;Sr2;Cu04 would be located in the CLP region. 

We conclude that the absence of a c-axis plasmon in the 
normal state of high- Tc cuprate superconductors results 
from an anomalously strong damping of the transport 
perpendicular to the Cu02 planes. This assignment also 
alleviates the discrepancy by two orders of magnitude 
existing in the literature between the plasma frequency 
obtained from optical experiments and that from the 
single-electron inter layer- hopping matrix element. The 
appearance of a c-axis phason mode, as observed in the 
superconducting state, is a direct consequence of the 
anomalously strong damping of the charge-carrier trans- 
port perpendicular to the planes. The decrease in a for 
T < Tc can be fully accounted for by the opening of a 
gap, the presence of which is a necessary consequence 
of the /-sum rule. There is no indication for a sudden 
decrease of the optical decay rate for E/ /c when T < Tc. 
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FIG. 1. Upper panel: Reflectivity of Lai.g5Sro.i5Cu04 
at 10 (squares), 20 (circles), and 27 K (triangles) with 
least-square fitted curves. Dashed curve: flt-parameters for 
8 K taken from Tamasaku et al. [Q]. Inset: Reflectivity 
at 33 K (triangles) with Drude-Lorentz fit. The parame- 
ters are e^c = 4.75, ctq = 7 f2~^cm~^, 7 > 4600 cm~^, so 
that Vp > 1200 cm~^. The TO phonons were parametrized 
with {uji (cm-^), 7i (cm-^). Si} (i = 1---4): {241.5, 24.8, 
17.6}, {312.0, 3.0, 0.019}, {351.0, 3.0, 0.044}, {494.1, 12.8, 
0.33}. Lower panel: Optical conductivity obtained from 
Kramers-Kronig analysis at 10, 20, and 27 K (solid symbols, 
from bottom to top). Chained curve: -Ime"^ X 13.00"^ cm 
Open symbols: conductivities with the chained curves sub- 
tracted. Dashed curve: calculation based on the Mat- 
tis-Bardeen model, taking into account the phonon contri- 
bution. 

FIG. 2. Temperature dependence of ctq (solid symbols) 
and (open symbols) corresponding to the fits in Fig. 1. 
Squares: sample 1. Circles and triangles : sample 2. Open 
and closed squares and open circles: with a constant ioo. 
Open and closed triangles: with full frequency dependence 
of phonon oscillators in ioo . Solid circles: conductivities at 50 
cm~^, after subtracting the loss function peaks as indicated 
in the lower panel of Fig. 1. Dashed curve: 1 — (T/Tc)*. 
Chained and solid curves: simulation for s- and d-wave pair- 
ing, respectively. Dotted curve: s-wave pairing plus constant 
background conductivity. The simulations of as/cr„ are at 
1/ = 50 cm^^. 



FIG. 3. Phase diagram in the {') / Vp)-{Eg / hvp) plane [ ^l| . 
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